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Abstract The momentum operator p = — i'V has radial component p = — ir(^ d r r) . 
We show that p is the space part of a 4-vector operator, the zero component of which is a 
positive operator. Their eigenfunctions are localized along an axis through the origin. The 
solutions of the evolution equation i dt ip = P° ip are waves along the propagation axis. Lorentz 
transformations of these waves yield the aberration and Doppler shift. We briefly consider 
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spin-half and spin-one representations. 



I. Introduction 

Our everyday experience of photons is that they are localized along their propagation axis, 



CO 

> 

s * ■ whereas in the usual formalism a photon of definite momentum (an eigenvalue of p ) is a plane 

wave spread over all space. The operator p = — i V can be split into a radial component p 
and an angular (transverse) component due to the identity 

O- i V = - i f (- d r r) - - - (f x L - L x f ) = p + p (1) 
r 2 r 

' with f = r/r and L the angular momentum operator. We note that 

Pk r 

pip(r) = -i-d r (rip(r)) (2) 
is undefined at the origin. In the appendix we show that continuity in r on opposite sides 
qh of the origin, i.e. that 



r ip(er) — r ip(— er) — > as e^O (3) 

is a sufficient condition for p to be a symmetric operator with respect to the usual inner product 
space 

(Vi(r)|^ 2 (r)}= /"d 3 rVi*(r)V 2 (r) = ( d£l I dr r 2 ^(r, 5,^)^(^5,0) . (4) 



a Jo 

The eigenfunction of p with eigenvalue k is 

Uk = — [<S(f,k) exp(iA;r) + S{— f, k) exp(— ikr)] (5a) 

' [<5(f,k) + *(-f,k)] exp(ik-r) (56) 



V 27r kr 

where k = k/|k| , k = |k| and 6(r, fc) is the delta function on the unit sphere (for this no- 
tation see the appendix). The equivalence between (5a) and (5b) may be seen by writing 
<5(?,k) exp(ifcr) = S(r, k) exp(zA;(f • f)r) = 8(r, k) exp(ik(k ■ f)r) = 8(r, k) exp(ik • r) , and 
similarly for 6(— r, k) exp(— ikr) . Then from (5a) 
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piik = — [rk5(r, k) exp(ifcr) — rk5(— r, k) exp(— ifcr)] 

[kfc<5(r, k) exp(zA;r) + kfc^(— f,k) exp(— zfcr) ] = kuk- (6) 



V27r fcr 

The state Uk is localized along the entire k axis r — Ak (—00 < A < 00), and the density 
(u k Uk) is evenly distributed along the k axis as 

/ ( u k u k) d s r = (b-a) /irk 2 . 

J r=a 

Note that 

Vu u =K = u {k ^-k} 

where V is the parity operator. 

The Uk satisfy the orthogonality and completeness relations 

J d 3 ru* k (r)u k ,(r) = 6(k-k'), J d 3 k u* k (r) u k (r') = 6(r - r') (7) 

which are verified in the appendix. As the Lorentz invariant measure is d 3 k/k rather than 
d 3 k , this suggests defining the inner product spaces 

d 3 k 



and the basis state 



fc(k)|&(k)) i/fc = J ^#(k)0 2 (k), (8) 
(^(r)|^(r)) i/p = J ^^(r)^(r), (9) 



Wk = Vkru-k = —j= [S(r, k) exp(ifcr) + 5(— f , k) exp(— ikr)] (10a) 
V27rfcr 

[ <*(?, k) + <&(-?, k) ] exp(i k • r) , (106) 



^/2nkr 

which is an eigenfunction of the operator 

p = - i r(-^= d r y/r) =-i r(d r + ^-) (11) 
^/r 2r 

which is symmetric with respect to the 1/r inner product space (9). From now on we will be 

concerned with the operator p and its eigenfunctions Wk (rather than with p, ) because tWk 

transforms as a scalar under Lorentz transformations (see below Sec. 3). From (7) 

<ufc(r)|uv(r)) 1/r = k6(k-k'), < Wk (r)| Wk (r')> 1/fe = rJ(r-r'). (12) 

Given a ^(r) which is a superposition of various % 

^ (rH / ^( k W=(<l<K k )) 1/fc (13) 
then the orthogonality relation (12) enables one to project out the distribution 4>(k) as 

S I <Kr) ) i/r = <V| I ^ 0(k') «v ) i/r = m ■ 
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Expanding out the quantity (w^ \i/)(r) ) we show that this is a Fourier transform of ip(r) 

\ /l/r 

as well as of r) along the line r = Ak, (—00 < A < 00), as 

w k I V(r)\ = 0(k) ee / ^V(r)<(r) 
/l/r J r 

— f drr J dQ, ip{ri) [#(r, k) exp(— zfcr) + <5(— f,k) exp(ifcr)] 

Jo Jfi V2nkr 

= ^_ y cZr y — [^(rk) exp(— ifcr) + rk) exp(ifcr)] . 

It is more transparent here to write k in spherical coordinates: with k — (k,$k,<f>k) then the 
above is 

4>{k,e k ,4>k) = y ^ ^ [V^A.^k) exp(-zfcr) + ip(r,n- 9 k ,<t> k + ir) exp(ifcr)] 
[(F c -if s ) + {F c + iF 8 )V]V^i>{r,9 k Ak) (14) 
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2Vk 
= Uip{r, 6 k ,(p k ) 



T c f(r, 6, 0) = g e (k, 6, 4>)= \l- I f(r, 0, 0) cos(rfc) dr , 



where V is the parity operator and J- c , J- s are the Fourier cosine, sine transforms dehned by 

T a f(r, 6, 4>) = g s (k, 6, <j>) = ^ J f(r, 6, cf>) sin(rfc) dr . 
The operator of (14) 

U= —7= \(f c - iF a ) + [T c + iT a )V\ VF (16) 
2\/k 

defines a unitary mapping from ip(r) to 0(k) whose inverse is 

U = -^[(T c + if s ) + (f c -if s )V]Vk (17) 

where the T c T c are as in (15) but with r, k interchanged. It is straightforward to check that 
UU — 1 , with the aid of T C T C — T s T s — VV = 1 . We emphasize that U, U are in a sense 
one dimensional transform operators: U maps ip(r) on any axis through the origin onto 0(k) 
on the same axis. 

As pu>k — ktUk , and k is the space part of the 4-vector k x = (k,k) , we look for the 
operator p° such that |5° w k = kw^ . (One can see by inspection of (10a) that % is not an 
eigenfunction of the "radial momentum operator" p r = —ir~ x l 2 d r r x l 2 which anyway is not a 
Hermitian operator [1]). The operator p° is multiplication by k in momentum space, i.e. 

p° =UkU (18) 

so that p° is a positive operator. To simplify (18) we need the following identities [2] 

J~ C J~ c — J~ S J~ s — 1 , J~ S J~ 2 — 7Ye , *F cJ~, s — 7~to (1^) 



where H e , H are the Hilbert transforms of even, odd functions defined by 

/(*Atf). + , 2 [°° /(Ar) 



d\. 



A 2 
'( 

A 2 



(20) 



We will write 

T± = T c ± i T s 

then from (19) there follows the further identities 

T+ — — i (H e - H D ) , F- F- = i (H e - H a ) , 

T+ = 2 - i (H e + H ) , T- T+ = 2 + i (H e + H ) , 

also 

T ± k = Tid r T ±7 kT ± = ±iT ± d r . (21) 

Returning to p° — UkU we have 

p° = UkU = ^= (T + + T-V) k (T_ + T + V) 

= i-^d r (-?+ + T-V) (T- + T+V) Vr 

= i-^d r + T-T+) + (T-T- - T + T + )V\ yft 

1 :8 r [ (H e + Ho) + (He-H )V\ y/¥ = - ^= d r H+^ (22) 



where 



V— I v *-e i ' i ' \ 6 ' ■'O ) 'iv' j— ' 



H+/(r) = ^ [ (H e + H„) + (W e -H )V\ /(r) (23a) 



- 1 /""//(Ar) _ /(-Ar) \ 
7T i V A - 1 1 + A / 
1 p /(Ar 



7T 



A-l 



dA 

dX (236) 



which is the Hilbert transform of f(r) along the axis Ar, (— oo < A < oo) . 

We now establish that jf = — drH+^/rw^ = kw^. For simplicity we will choose 
the particular case when k = (0, 0, k) , k > . Then 

w k 3 = 1 [ 5(f 3 - 1) + 5(-f 3 - 1) ] exp(i k z) 
\j2ixkr 



and 



H + exp(i k z) = +i exp(i z) z > 

— i exp(i fc z) z < 

-^=H+Vrw k 3 = * = [5(f 3 - 1) - 5(-f 3 - 1)] exp(ifcz) 

-^d r H + y^w k 3 = 1 fcf 3 [o(f 3 - 1) - o(-f 3 - 1)] exp(ifcz) 
V r V2nkr 

: fc [ £(f 3 - 1) + <5(-f 3 - 1) ] exp(i kz) = k w k 3 



\f2nkr 

The operator components <9 r , H+ of p° do not commute. If in the working out of (22) we take 
k to the right instead of to the left, we arrive at 

p° = UkU = - ^j= [ (H e + H ) - (H e -Ho)V) d r ^ 

= -^=H-d r Vr (24) 



where 

H_/(r) = X - [ (H e + H ) - (H e - H ) V ] /(r) (25a) 
= - !°° S gn(X)^-dX. (256) 

As 

= Ti+H- = — 1 , 

which can be verified their definitions (23a), (25a) together with H e H — H Ti. e = — 1 , then 

P° 2 = 4= 9r H + H-d r = -4= % = P 2 

as necessary. Also with the adjoint relations (-^=H" y/r)^ = — -±=H° \/r then from (23a), 
(25a) 

2. Wave equations 

We can write a wave equation substituting p for the usual p : 

[(id t f - p 2 ]^ = o 

[_ a 2 +r -l/2 a 2 r l/2 ]v , = (2g) 

which is essentially the wave equation —\Jnf> = [—d 2 + \d 2 r — \ L 2 \ip — with the operator 
component \ 1? excluded. This modified wave equation (26) has solutions if,^, V'k which are 
eigenfunctions of p : 

■i/>k = exp(-ikt) tWk, — exp(ifci) w£ . (27) 

Both the states V'kiV'k are waves localized along the propagation axis r = Ak (— oo < A < oo) 
and proceeding in the +k direction, and accords with our everyday experience of photons being 
localized along their propagation axis. If if, satisfies (26) then there is a conserved density 

a = i(if)*d t if)-(d t r)i> 

However this density is indefinite. 

We will from now consider the first order evolution equation 

id t ifj = fif, = - r- 1/2 d r H+ r 1/2 if; . (28) 

Now only the positive energy eigenstates ■i/'k of (27) are solutions and the negative energy 
components if>^ are excluded. Consider the non-negative density p 

p=±[iP*if)+ {r- 1 ' 2 H + r 1 / 2 if)*) (r-^ 2 H+ r x ' 2 if,) ] 
r 

then 'dip is a solution of (28) 

d tP = - d t if)* if) + - (r- l ' 2 H+ r 1 ' 2 if,*) (r- l ' 2 H+ r 1 ' 2 d t if)) + CC 
r r 

= -- {r- l ' 2 d r H+ r 1 / 2 ip*)if) + - {r- x l 2 H+ r 1 ' 2 if)*) (r" 1 / 2 ^ H- d r r x l 2 if)) + CC 
= - i {r- 1 d r H + r 1 / 2 if)*) (r- 1 ' 2 if)) - i (r- l H+ r 1 / 2 if)*) (r^r 1 / 2 if,) + CC 



V f 



i(r- 1 W + r 1 /V)(r- 1 / 2 V) + CC 



V • J 
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where CC stands for the complex conjugate terms and we have used the operator identity 

(r- 1 d r rA)B + A{r- l d r rB) = V • [f AB}. 

Substituting ^ = into the above expressions for the density and current, both the density 
and current are evenly distributed along the k axis r = Ak (— 00 < A < 00) , and the current 
is uniformly in the + k direction. 

3. Lorentz transformations 

The momentum space transformation is straightforward. The fc A — (k,k) = (k, k 1 , k 2 , k 3 ) 

transforms in the usual way as a 4-vector, that is under a finite boost of velocity v in the z 
direction then 

k' 1 = k 1 , k' 2 = k 2 , k' 3 = 7 fc 3 - ~/vk , k' = 7 fc - -yv k 3 , (29) 

and 

k' 1 = jfeV7(l - vk s ) , k' 2 = fc 2 / 7 (l - vie 3 ) , k' 3 = (jfe 3 -v)/(l-vk 3 ) (30) 

with 7 = (1 — v 2 )^ 1 / 2 . The last of (30) is the well known aberration formula, as putting 
k 3 — — cos 6 (light observed from the direction k has momentum oc — k) then 

cos 9' = (cos 9 + v)/(l + v cos 9) . 

The basis state — exp(-ikt) % transforms as a scalar in momentum space, i.e. 

(W = (^k)k^k' = [S(-f, k>) + 6(r, k>) } expHfc'i + ik' ■ r) (31) 

V zirk r 

where k' is given by (29). From (31) can be inferred the frequency change from the exponential 
term, and also the change of the propagation direction k which agrees with the aberration 
formula. 

The k transformation is generated by the boost and rotation operators in momentum 
space which are 

N fe = i k (d/dk) , L fe = - i (k x d/dk) (32) 

then (29) is k' — exp(i> A^ 3 ) k . In configuration space the boost transformation will be compli- 
cated due to the fact that jf is a non-local operator, so there is no simple transformation of 
the configuration space r coordinates corresponding to (29). The Lorentz boost generator N 
must satisfy 

[N,p°]=ip [N\p b ]=i5 ab TP (33) 

which implies that N as well as jp is a non-local (integral) operator. The operators generating 
the Lorentz transformations in configuration space are 

N = (r- x l 2 H- r 1 ' 2 ) (r V - 2 f d r r) = (r V - 2 f 8 r r) (r- l ' 2 H+ r 1 ' 2 ) 

(34) 

L = — ir x V 
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which are Hermitian with respect to the 1/r inner product space (9). That the N defined 
above satisfy the relations (33), and that the Lorentz group relation [N a ,N b ] = —ie abc L c is 
also satisfied we leave to the appendix. An infinitessimal boost <5v causes a change in the 
wavefunction ip — > (1 + iSv ■ N)^. In the appendix we calculate Nwk and show that the 
combined transformation N u>k + Nfc Wk — , which means that is a Lorentz scalar. 

This completes our investigation of the spin-zero operators p° , p , N , L which together 
obey the Poincare group commutation relations. 

4. Spin-i and spin-1 representations 

Massless spin-i and spin-1 Hamiltonians can be constructed which are differential operators, in 
contrast to the non-local p° Hamiltonian for the spin-zero case. The spin-i massless (neutrino) 
Hamiltonian equation 

id t V = <T-piH = -i(<T-r)(d r + ^-)y (35) 

has as eigensolutions either column of the 2x2 matrix e ±lkt (er • k ) u>k • Formally one can 
also construct a spin-^ Dirac-like equation with mass idt^ = (a. ■ p + f3m) , which has as 
eigensolutions a column of the 4x4 matrix [exp(— i(m 2 + fc 2 ) 1 / 2 t)(a • k + (3 m) Wk] , however 
in the rest state when k = then is undefined. 

The source- free electromagnetic field is a massless spin-1 field. Defining the complex field 
F = E + iB , then the source-free Maxwell's equations are 

[dt + i V x ] F = , V • F = . 

which we will replace by 

[d t - p x ] F = , p • F = . (36) 

A solution of (36) is 

F = e- ikt (w k3 , i w k 3 , 0) 
which is a wave localized along the z axis proceeding in the +z direction. 

5. Outlook 

The attraction of the p x operators is that their eigenfunctions conform with our experience 
of photons: that they appear to be localized along their propagation axis. On the other hand 
the chosen origin of coordinates plays a more significant role in this formalism - indeed a 
shifted origin (or observer) may not register the localized wave at all. In the usual case the 
exponentiated momentum operator {exp(a • p) } has a clear geometrical role: in shifting the 
coordinates by a . By contrast the meaning of the operator {exp(a ■ p) } is obscure. 

Mathematically the the p x operators are interesting because they demand continuity 
through the origin, in contrast to the usual radial operators whose use is limited apart from 
describing sources or sinks at the origin. 
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Appendix 

That the operator p is symmetric with respect to the inner product (4) 
Expanding the inner product 

0(r)|p^(r))= ( dn [ dr r 2 0*(r) [- irr- 1 d r r^{r)} 
' Jn Jr^O 

r oo 

dr r 2 [— i fr _1 d r r<j)(r) ]* ip{r) —if dflr (rcf)*(r)) (rip(r)) 

Jn 

= (p0(r) U(r) ) - i [ dftr [(r0*(r)) {rip{r))] (Al) 

assuming the (r</>), (rip) tend to zero sufficiently fast at infinity. But with the condition (3) 

r ip(er) — r tp(— er) — > as e^O 

then the last term on the RHS of (Al) is zero, as the contributions form opposite sides of 
the origin cancel. Thus p is symmetric. Note that it is the presence of the f within the 
p operator that leads to symmetry, by contrast the operator p r = —ir~ l d r r has the more 
onerous requirement for symmetry that rip(r) ^ as r — > 0, which is the underlying 
reason why p T is not Hermitian (as discussed by Hell wig pi 76 [3]). 

The delta junction on the unit sphere 
In spherical coordinates 

f = (sin 9 cos <f> , sin 9 sin </> , cos 9 ) 

k = (sin 9k cos 4>k , sin 9k sin <f> k , cos 9k ) 

then 

5(r , k) = <5(cos6> - cos0 fe ) 5(cp - <p k ) (A2) 
8(-r , k ) = <5(cos 9 + cos 9 k ) 5((f> + tt - cj> k ) . (A3) 

Orthogonality and completeness of the 
Expanding out (7a) 

/d 3 r u£(r) Uk'(r) — [ dQ I dr — \— [S(— r,k) exp(ifcr) + S(r,k) exp(-ikr)] x 
Jn Jo 2nkk' 

[6(— r, k') exp(— ik'r) + <5(r, k') exp(ifc'r)] 

= -i- <&(£, k') 5(jfe - fc') = 5(k - k') . (A4) 

The boost operator N of (34) 
First consider the operator 

N'E-j(rV- 2f <9 r r) (A5) 

and defining 
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the following commutation relations may be directly verified 

[N 7 ,a°]=is, [N' a , s h ]=i5 ab s° . 



is, [iv~ , s" j — io~~ s" . (A7) 

Also the operator —id r r commutes with N' so that 

= -ir-^=d r y/?) = (i-^=d ry /?, p) (A8) 

V r v r v r 

also obeys (A7) with t x instead of s x . Finally we multiply t° on the right by (i-j=H + \/r) and 
N' on the left by (i-^=H- y/r) to make the new operators p° , N , which recalling H+H- — — 1 
also satisfy the relations corresponding to (A7), i.e. 



N = (lH-V?)(rV - 2?d r r) =(rV - 2r d r r) (-^H+Vr~) 



-4= d r H + \fr — ^= W_<9 r v^, 



-ir—=d r \fr . 



Finally the fact that 



[N a , iV 6 



-t e abc L c 



(A9) 



(AW) 



follows from [N' a , N' b ] = -i e abc L c where L = - i r x V . 



To s/iow i/mi NiOk + Wk = 

We will consider the boost operators A 3 , N% of (32), (34), and recalling from (10) we have 
in spherical coordinates 

A| — i cos 6 k k d k - i sin 6 k dg k 

N 3 = (- cos6rd r - sm6d e ) (^=H + y/r) 



Then 



and 



1 



w k = 



2\lixkr 



[8(cos6 



h) 6~{(t> - 4>k)) ex.p(ikr) 



+ 5(cos 6 + cos ) 5(0 + 7r — 0^ ) exp(— ikr) ] 

1 i 

(—=H+\/r )wk = — [<5(cos6> - cos^ ) 8(<j> - 4>k)) exp(ikr) 

V r 2\/7ikr 

— 5(cos 6 + cos 6 k ) 5(4> + n — cf>k) exp(— ikr) ] . 

exp(ikr) 



Nj: «k = + i cos 6k 8(cos 6 — cos 6k ) 8(<p — <pk ) k dk [ 



2v irkr 



- i t= 8(4> - 4> k ) ) sintffc dg k [<5(cos6> - cos 6^ )] 



2V7rfcr 

+ i cos 6k 8(cos 6 + cos 6 k ) 6(<fr + it — 4> k ) k d k 



|.exp(— ikr) 



- i o~(<t> + tt - (f> k ) ) sin^fe dg k [8 (cos 6 + cos6 k )} 



2-Jirkr 



N 3 Wk — — i cos 6 8(cos 6 — cos 6 k ) 8(4> — (f>k)r d r [ 



exp(ifcr) - 
2\fnkr 

- i exp f^ S((f> -<f> k )) sm6d [S(cos6 - cos6 k )} 



2Vnkr 

+ i cos 6 8 (cos 6 + cos 6 k ) 8(4> + tt — 4> k ) r d r [ 



exp(—ikr) 
2\fr7kr 



+ i exp ^_^_ r ) s^ + n-^k)) sin 6 dg [8(cos 6 + cos 6 k )} . 



2\/iTkr 

Now if we add (A12) and (A12) together then the components cancel in order. 



(All) 



(A12) 
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